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We show that the wave packet of a biphoton generated via spontaneous parametric down conver- 
sion is strongly anisotropic. Its anisotropic features manifest themselves very clearly in comparison 
of measurements performed in two different schemes: when the detector scanning plane is perpen- 
dicular or parallel to the plane containing the crystal optical axis and the laser axis. The first 
of these two schemes is traditional whereas the second one gives rise to such unexpected new re- 
sults as anomalously strong narrowing of the biphoton wave packet measured in the coincidence 
scheme and very high degree of entanglement. The results are predicted theoretically and confirmed 
experimentally. 
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INTRODUCTION 

Quantum entanglement is one of surprising conse- 
quences of quantum mechanics. It lies in the center of 
attention since the famous paper by Einstein, Podolsky, 
and Rosen[lJ (EPR). Two (or more) subsystems are en- 
tangled if the system as a whole is characterized by a 
wave function or a density matrix, which cannot be pre- 
sented in the form of products of subsystems' wave func- 
tions or density matrices Q. In some cases entangle- 
ment (for example Bell states) leads to complete dctcr- 
minacy of the combined bipartite system while states of 
the subsystems are fully undetermined. Entanglement is 
the key element of such phenomena and domains of mod- 
ern quantum optics as quantum teleportation, quantum 
cryptography, Bell violation experiments, quantum com- 
putation. 

A special class of entangled systems, attracting a per- 
manently growing attention, is that of systems with con- 
tinuous variables. The most often and widely studied 
example is Spontaneous Parametric Down Conversion 
(SPDC) 0, [JS S, 0, !, 0, E^. Some questions arising 
in discussions of entanglement in such systems are how 
high is the entanglement in a given multipartite state, 
what is the best measure of the degree of entanglement, 
how the degree of entanglement can be measured exper- 
imentally, etc., etc. It's rather well recognized [1, [ll| 
that in the case of pure bipartite states the degree of 
their entanglement can be characterized by the so called 
Schmidt number K (T^. [l3| . However, the Schmidt num- 
ber definition does not provide any recipes on how it 
could be measured. An alternative approach 1(J 14 1 to 
the entanglement characterization is based on analysis 



of conditional and unconditional probability density dis- 
tributions, which are related directly to experimentally 
measurable coincidence and single-particle wave pack- 
cts of particles. Specifically, it was suggested [14J to 
characterize the degree of entanglement by the ratios 
R x = Ax^ I Ax^ or R k = Ak^/Ak^ c \ where Aa:0' c ) 
and Afc( s - c ) are the single-particle (s) and coincidence (c) 
coordinate and momentum wave packet widths, x and k 
are the coordinate and momentum variables of one of 
two particles under consideration. As it was shown [Icj |. 
K = R x = Rk for a rater general class of double-Gaussian 
bipartite wave functions 
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where x\ and x 2 are continuous variables of two parti- 
cles, a, (3, 7, S and a, b are constants. Physical meaning 
of these parameters is determined by a problem under 
consideration. There are problems in which a double- 
Gaussian bipartite wave function arises quite naturally. 
One example of such a phenomenon is the strong-field 
electron-positron pair production [l5j |. In some other 
processes and phenomena the arising bipartite wave func- 
tions arc not double-Gaussian but can be successfully 
modeled by double-Gaussian functions with appropri- 
ately chosen parameters. In these cases the parameters 
K, R x , and Rk are usually close if not equal to each 
other. And, finally, if a bipartite wave function has a 
form strongly different from the double-Gaussian one, all 
three parameters K, R Xl and Rk can differ from each 
other. We assume that in such a case each of these pa- 
rameters characterizes the degree of entanglement to be 
seen in different ways, and entanglement itself cannot be 
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characterized completely by any single universal param- 
eter. 

Note that for non-entangled states K = Rx = Rk = 1. 
Moreover, for non-entangled states coincidence and single 
wave packets are identical to each other. Any differences 
between the coincidence and single wave packets can be 
considered as a manifestation of entanglement. Macro- 
scopic characteristics of these differences are given just 
by the parameters R x , Rk, and K. If these parameters 
differ from each other, the largest of them indicates an 
optimal way of observing entanglement experimentally 
or calculating mathematically, whereas in other ways of 
measurements or calculations entanglement can be hid- 
den. 

Also, one can characterize the degree of entangle- 
ment by the EPR-rclated parameter given by the in- 
verse coordinate and momentum conditional uncertain- 
ties or wave packet widths for any given particle: Cepr = 
l/[Ax^ x AfcM] [lEEEl]. For non-entangled states a 
specific value of this parameter depends on the definition 
of widths and on the shape of single-particle wave func- 
tions ip\{x\) and ^(a^), the product of which gives the 
non-entangled bipartite wave function. If the widths are 
determined as square roots of variances (as in the work 
Q), and the functions ipi(x\) and "02 (^2) are Gaussian, 
the no- entanglement value of the EPR-relatcd parame- 
ter is given by C E p R var = 2. If the widths Ak^ and 
Ax^ are determined as the widths of the correspond- 
ing probability distribution curves at the half-maximum 
level (as in this work), for the same forms of 1^1 and 
V> 2 , Cepr 1/2 = (41n2)- x w 0.36. In terms of the 
EPR-rclated parameter Cepr, the degree of entangle- 
ment for non-factorized bipartite functions is evaluated 
as Cepr/C e pr . 

It is assumed usually that deviations of realistic bipar- 
tite wave functions from the double-Gaussian one [Eq. 
(fTj)] do not give rise to any pronounced new features 
of the processes under consideration and, in particular, 
all the entanglement parameters remain of the same or- 
der of magnitude if not equal identically. In contrast 
to this assumption, we show below that, in the case of 
SPDC, deviations of the biphoton wave function (to be 
found below) from the double-Gaussian one are well pro- 
nounced and they can change drastically the observable 
wave packet pictures. Owing to this, the SPDC entan- 
glement phenomenon appears to be essentially multipara- 
metric, with the above-indicated parameters significantly 
differing from each other. 

Concerning the observable biphoton wave packet pic- 
tures, formulated explicitly or not, the usual notion is 
that they are more or less invariant with respect to rota- 
tions around the pump laser axis. In contrast to this, we 
show that the SPDC biphotons are strongly anisotropic 
and the degree of anisotropy depends essentially on rel- 
ative orientation of a crystal and the observation di- 



rection. Furthermore, usually the biphoton momentum 
wave packet structure is supposed to be controlled by 
the laser pump angular distribution and by geometrical 
factors like crystal shape etc. 0, S S]- That is why it 
seems that in order to generate well-localized biphoton 
wave packets in the fc-domain one has to use good laser 
sources with a very small angular divergency and/or long 
crystals. At the same time forming particular angular 
profile of the pump is needed for coupling biphotons into 
the optical fibers 3, 17 1, that is very important for some 
applications. So the question arises whether it's possible 
to get a good biphoton wave packet localization with a 
poor laser having a relatively high angular divergence? Is 
it possible to vary the degree of entanglement while keep- 
ing constant parameters of the pump and nonlinear crys- 
tal generating photon pairs? The answer we give in this 
paper is yes, owing to the above-described anisotropy, 
it's possible to control the degree of wave packet local- 
ization and the degree of entanglement simply by tuning 
a crystal orientation with respect to detectors position. 
Under the optimal conditions the achievable degrees of 
wave packet localization and entanglement appear to be 
very high. In particular, it becomes possible to get the 
coincidence angular biphoton distribution much narrower 
than that of the pump, which can look rather unexpected 
and unusual from the point of view of traditional theoret- 
ical and experimental approaches to the investigation of 
SPDC entanglement. The results outlined above follow 
from the theoretical formulas to be derived below and 
they are confirmed by direct experimental observations. 



THEORETICAL DESCRIPTION 
Derivation of the main general formula 

Let us consider a collinear and degenerate type I SPDC 
process, when extraordinary pump photon of a frequency 
Up decays in two ordinary photons (signal and idler) 
with equal frequencies lu p /2 and propagating more or less 
along the pump beam. In 3D, refractive index surfaces 
for extraordinary [n e (r; oj p )} and ordinary [n (f; u p /2)] 
waves in an anisotropic crystal are, correspondingly, an 
ellipsoid and a sphere. Fig. [TJa) shows sections of these 
figures by three coordinate planes (xz), (xy), and {yz). 
The optical axis of a crystal is taken directed along the 
Ox axis, and the orthogonal axes Oy and Oz are cho- 
sen in such a way that the (xz) plane contains the pump 
laser axis Oz'. The arc DB in Fig. |Tfa) is a part of the 
circle by which the sphere and ellipsoid cross each other. 
We assume that the laser axis is directed strictly to the 
point B where the arc DB crosses the xz plane. This 
means that for pump and emitted photons propagating 
strictly along the z'-axis the collinear degenerate phase 
matching condition is exactly satisfied, k p = k\ + lt2 at 
k p || fci || ki || Oz' , ki = &2 = kp/2. In our consideration 
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FIG. 1: (a) Octant of the refractive index surfaces n e (f) and n (r) 
for pump and signal photons, (b, c) are two perpendicular section 
by the planes _L and || shaded in (a); DB^ is the projection on the 
_L-plane of the circle shown partially by the arc DB in (a); i? p is 
the angle between k p and the z 1 axis. 



we assume that the pump is not a single plane wave but 
is given by a coherent superposition of plane waves with 
wave vectors k p filling in a cone, the axis of which co- 
incides with the Oz' axis, and the angular width a is 
finite. 

0' in Fig. [2(a) is some point in a far zone located 
at the laser axis Oz' and such that detectors registering 
photons are located in the plane perpendicular to Oz' 
and containing the point O' . The axes O'x' and O'y' in 
this plane are perpendicular to each other and to O'z', 
and O'x' G (xz), 0'y'\\Oy. Let us assume that detectors 
are installed along some line 0'£ in the plane (x'y') and, 
hence, they register only photons (1) and (2) with wave 
vectors k\ and k 2 belonging to the plane (z'£). Orien- 
tation of the detector-installation line 0'£ in the plane 
(x'y') can be changed with two limiting cases given by 
0'£ = 0'y'\\Oy and O't = O'x' € (xz) . In these two lim- 
iting cases the detectors register only photons with wave 
vector fci and k 2 belonging, correspondingly, to the yOz' 
and xOz planes. In Fig. [TJa) these planes are shaded and 
they are labeled with the symbols _L and ||, respectively. 
Here "perpendicular" and "parallel" mean that the ob- 
servation plane (i.e., the plane containing wave vectors 
of photons to be observed) is perpendicular or parallel 
to the plane containing the crystal optical axis and the 
laser axis. Below these two cases are referred to as those 
of the perpendicular and parallel geometry. Sections of 
the 3D refractive index surfaces by the _L and || planes 
are shown in the pictures (6) and (c) of Fig. [T] 

Note that the idea of changing orientation of the 
detector-installation line is used here for simplification 
of illustrations like those given in Fig. [TJ In real ex- 
periment the laser pump axis was horizontal (|| Oz), and 
detectors were installed in such a way that they registered 
only photons with horizontally oriented wave vectors k\ 
and k 2 . This was a crystal that was rotated around the 
laser axis (together with the laser polarization) instead 
of changing detector positions. There are the following 
rules of one-to-one correspondence between these experi- 
mental conditions and the above-described geometries of 
Fig- El (1) orientation of the crystal optical axis in the 
horizontal plane (yz) is equivalent to the || geometry and 



(2) orientation of the crystal optical axis in the vertical 
plane (xz) is equivalent to the _L geometry of Fig. [T] 

The second note concerns an additional assumption to 
be used below for simplicity and not motivated by orien- 
tation reasons. In terms of notations used in Fig. [lja), 
let us assume that the photons (1) and (2) with wave 
vectors k\ and k 2 belonging to the plane (z'£) arise only 
owing to the decay of pump photons with wave vectors 
k p belonging to the same plane (z'£). Rigorously, this 
is not necessarily true: owing to the non-collinear phase 
matching processes photons (1) or (2) can be emitted 
with wave vectors k\ or k 2 belonging to the plane (z'£) 
even if the pump wave vector k p is oriented differently. 
Such processes can affect the single-particle momentum 
distributions of emitted photons, and this effect is dis- 
cussed below. But, on the other hand, the non-collinear 
phase matching processes can be effectively suppressed 
by means of axially asymmetric focusing of a laser beam. 
This possibility is also discussed below and its experi- 
mental realization is presented. 

Under the formulated assumptions, with the help of 
the results derived by Monkcn ct al in 1998 [f|, in the 
approximation of a wide crystal we can write down im- 
mediately the following expression for the biphoton wave 
function depending on the transverse components of the 
emitted photon momenta k\ £ and k 2 £ 
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where sinc(u) = sin(u)/u, L is the crystal length in the 
z'-direction, and A z > is the longitudinal detuning 
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The above-mentioned approximation of a wide crys- 
tal means that its size in the direction parallel to the 
0'£ axis is large enough to make the transverse mo- 
mentum conservation rule, or transverse phase match- 
ing condition, exactly satisfied to give k p ^ = ki ^ + k 2 £. 
This relation is used in the definition of the momentum- 
dependent pump field-strength amplitude in Eq. 
E*(kp£) = E*(k\£ + k 2 ^), and it can be used also in 
the definition of the longitudinal detuning A z / ([3]). 

To simplify further Eq. ([3]) for the longitudinal de- 
tuning we can use the near-axis approximation in which 
1^1,2 {| *C fci,2 f k P £ <C k p . By expanding square roots in 
Eq. ([31) in powers of transverse components of all the 
wave vectors and keeping only two first orders we get 



A z ' = k p — k\ — k 2 



(4) 



where in the second term on the right-hand side we took 
k p w 2k\ = 2k 2 . As for the first term, k p — k\ — k 2l if it can 
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be taken equal zero identically, then Eqs. |[5J) and (|4|) give 
immediately the well-known and widely used formula of 
the work [6| 



)*E* p (k 



S2U 



M 2 



4fc„ 



.(5) 



This is the formula that is most often met and used for 
theoretical description of the SPDC biphoton entangle- 
ment. Eq. ([5]) has a form very close to that of the double- 
Gaussian wave function given by Eq. |T]). Indeed, very 
often the pump envelope has a Gaussian form. As for 
the sinc 2 -function, it can be rather successfully modeled 
by a Gaussian function too [1] . The key point of such a 
reduction to the double Gaussian form is the separation 
of variables k\ j and k 2 £ for their two independent linear 
combinations k\ £ + k 2 £ and k\ £ — k 2 ^ in the arguments 
of two protofunctions of Eq. (J5]), E* and sine. As we 
will show now, such separation of variables and reduc- 
tion of the wave function to the double-Gaussian form 
essentially depend on the assumption k p — k\ — k 2 = 
which is true only under very special conditions. 

As it can be seen from the pictures of Fig. [l] the 
collincar phase matching condition k p — k\ — k 2 = can 
be taken satisfied for different orientations of the pump 
wave vector k p only if the detection direction C£ is per- 
pendicular to the plane containing the laser and crystal 
optical axes, i.e., if 0'£\\Oy. In this case, for all wave 
vectors k p located in the _L-plane of Fig. QJa) and not 
too far deviating from the Oz' axis, their length can be 
taken approximately equal: kp^p) « fc p (0) = k\ + k2, 
where $ p is the angle between k p and Oz', |$| <C 1. Ow- 
ing to the last condition the arc DB± in Fig. can 
be approximated by an arc of a circle with the center at 
the origin O, and then equality of all wave vector lengths 
k p (i!}) becomes evident. 

At all other orientations of the detector-installation 
line 0'£, non-parallel to Oy direction, the difference 
kp{dp) — ki — k 2 on the right-hand side of Eq. ([4]) can- 
not be taken equal zero at $ p even approximately. 
At small values of i? p the dependence k p (i) p ) can be lin- 
earized to give 
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where n' p (ip) = dn p ((p)/dip, ip is the angle between k p 
and the crystal optical axis (Ox in Fig. [T|) and Lp is the 
angle between the optical and laser axes. On the other 
hand, as the laser axis Oz 1 and the detector-installation 
line 0'£ are perpendicular to each other, at small $ p we 
have k p (0) "d p k p £ = fci £ + k 2 £ . By remembering that 
kpifip = 0) = k\ + k-i and combining all these relations 
together we get the following expression for the longitu- 
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and for the momentum biphoton wave function 
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where n' p = n' p (ip ) and n Q Tip (ipo). Eq. © is the 
main new formula for the biphoton wave function we 
suggest instead of the traditional one given by Eq. ((5]). 
These two formulas differ from each other by the first 
term in the sum in square brackets under the symbol 
of the sine— function in Eq. ([5]). Below we will analyze 
how important this addition is. But before doing this, 
let us rewrite Eq. ((8]) in terms of the scattering angles d\ 
and 02, defined outside the crystal as 6>i i2 = 2fci :2 z/k P °\ 
where k p ^ = uj p /c. By using these definitions we express 
the transverse components of the emitted photon wave 
vectors via the scattering angles and reduce Eq. ((8]) to 
the form 
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where E p (9 p ) is the pump amplitude angular distribution 
outside the crystal, p = \{9\ + 02 ) = n a d p (to remind, 
by definition $p is the angle between k p and the laser 
axis Oz' in a crystal) . In terms of scattering angles the 
longitudinal detuning of Eq. ([7]) appears to be given by 
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The new first term in the square brackets of Eqs. © 
and (flTj]) is linear whereas the second term is quadratic 
in small angles 0i and 02- For this reason we can expect 
that the linear term is even more important than the 
quadratic one, if only the refractive index derivative n' p 
is not negligibly small. But this docs not mean that the 
quadratic term can be dropped because in such an ap- 
proximation we would get infinitely wide single-particle 
distributions. So, both linear and quadratic terms have 
to be taken into account. The role of the linear term 
is illustrated by three curves of Fig. [5] where the func- 
tion sinc 2 (LA z //2) is plotted in its dependence on 0\ at 
02 = and at three different values of the refractive in- 
dex derivative. The curves of Fig. [5] are calculated for 
the pump and crystal parameters corresponding to the 
experimental ones (see the following Section): the pump 
wavelength X p = 325 nm and LUO3 crystal of a length 
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sine 2 (0i, 6 2 = 0) 




0.02 0.02 0.04 



(c) 



0.01 0.565 0.575 



FIG. 2: The function sinc 2 (LA z //2) at 2 = and (a) n' p 
0, (b) - 0.01436, and (c) - 0.1436; 0i is in radians. 



of the sine 2 function does not give contributions to the 
single-particle photon momentum distributions (see the 
derivation and discussion below). The large-angle second 
peak of the sinc 2 -function was not observed also in the 
coincidence measurements described in the following Sec- 
tion simply because this was out of the detector scanning 
range. For these reasons we restrict our further analysis 
by considering only one peak of the sinc 2 -function located 
at small values of the scattering angle Q\ in both cases of 
zero and nonzero values of the refractive index derivative 



L = 1.5 cm. In this case we have found ipo ~ 60° and 
n' v = —0.1436 for the detection of photons in the vertical 
plane of Fig. [TJ The three curves of Fig. [2] correspond 
to (a) n' = - photon propagation and detection in the 
horizontal plane, (b) n' p = —0.01436 - some intermedi- 
ate case between _!_- and | -geometries of Fig. QJa), and 
(c) n' — —0.1436 - propagation and detection in the ||- 
plane. One can see, that with a growing value of \n' p \ 
the structure of the curves in Fig. [2] changes drastically. 
A single wide peak splits for two peaks, spacing between 
them grows, and the peaks themselves are getting very 
narrow. The widths of the only peak at n' p = and of 
narrow peaks at n' p = —0.1436 are equal to 24 mrad and 
0.5 mrad, correspondingly. The ratio of these two widths 
equals to 48 and this number characterizes the degree 
of peak narrowing arising owing to modifications in the 
formula of Eqs. ©, © compared to that of Eq. (f5|). 

Mathematically splitting of a single peak for two ones 
follows from the quadratic dependence of the detuning 
A z < (|10p on 6\. The sine— function is maximal at zero 
value of its argument, i.e., at zero detuning A z /. In a 
general case the quadratic equation A z /(0i) = has two 
solutions which correspond to two peaks in the depen- 
dence of sine 2 on Q\ . The only exception is the degenerate 
case n' p = when two solutions of the quadratic equa- 
tion merge into one. Qualitatively the same conclusions 
and the same arguments are illustrated by the Fig. W[c). 
In this picture the vector ki is plotted along the optical 
axis Oz' (0 2 = 0), and it ends at the point A. The end- 
ing locus for vectors k\ is given in this case by a circle of 
the radius no/2 (in units of oj p /c) and with the center at 
the point A. As it's seen from the picture, there are two 
points B an C where this circle crosses the ellipse n p (9 p ) 
which, in its turn, is the ending locus for the pump wave 
vectors k p . This picture shows that in the case n' p ^ 
there are two directions of the vectors k\ and k p (at a 
given ^2) in which the exact phase matching condition 
k p = ki + &2 appears to be fulfilled. Quantitatively, the 
nonzero angle of the second exact phase matching direc- 
tion appears to be very large. The value of 9\ ps 0.57 in 
Fig- El^c) corresponds to about 30° for the scattering an- 
gle 6*i and to 9 p f=a 15° for the pump. As the last value is 
much larger than the pump angular divergence (typically 
a ~ 4 mrad f» 0.23°), the second (nonzero-angle) peak 



The quadratic dependence of A z / on 6\ and a transition 
from the case when the equation A z ' = has a single 
solution to the case of two different solutions (at n p ^ 0) 
explain also qualitatively the reasons of peak narrowing. 
Roughly the sine 2 peak widths can be evaluated from 
the condition |A Z /| ~ 1/L. In the case n' = at 02 = 



we have A z i ~ k p ^9\ which gives A#i ~ 1/v/ Lk p °\ In 
contrast, in the case n' p ^ in a vicinity of each narrow 
peak the detuning is approximately linear in 6\ , and from 
the same condition |A 2 /| ~ 1/L we find A0i - l/(Lk p 0) ). 
As the product Lk p °^ is rather large (~ 10 5 ), both width 
are small but the peak widths occurring in the case n' p ^ 
is much smaller than that occurring in the case n' = 0, 

i/(Lk p o) ) « vv^fcF ■ 



Coincidence and single-particle distributions 

Coincidence distributions of photons are given sim- 
ply by the squared absolute value of the wave function 
^(#1 j 82) °f Eq. |9]) at a given value of one of the angles, 
0i or 2 . E.g., at 2 = 
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Single-particle distributions are given by the squared 
wave function of Eq. © integrated over, e.g., 02- In 
the case n' p ^ it is convenient to make the integration 
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of 0i ) to get 
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At ^ the sinc 2 -function is not negligibly small only 
in small vicinity of points where the detuning A 2 <(0; 0i) 
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turns zero. From the condition A z > = we get a 
quadratic equation in 9, solutions of which are given by 



(14) 



Only one of these two solutions (9 a ) is small enough for a 
region around this point to give a non-zero contribution 
to the integral over 9 in Eq. (|12p . Near the point 6 a the 
detuning A z i can be approximated by a linear function 
of 9 
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and the sine 2 function can be approximated by the delta 
function 
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to give finally the following simple expression for the 
single-particle momentum distribution of photons 



dw[ s) (9 1 ) 



d6x 



oc 



E,, 



,12 



2n'9i 



(17) 



with 9 a (9 1 ) given by Eq. (Tj 

The coincidence (solid lines) and single-particle 
(dashed lines) photon momentum probability distribu- 
tions are plotted in Fig. [3] for two cases, (a) n' p = (_L- 
geometry) and (b) n' p ^ (||-geometry). In both cases the 

dotted- line curves describe the pump intensity \E p (9 p )\ 2 
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FIG. 3: Coincidence (solid lines) and single-particle (dashed lines) 
photon distributions, dw^ /d9\ (at 82 = 0) and dw^ /ddi, calcu- 
lated for the cases (a) n' p = and n' p = —1.436. The dotted-line 
curves show the pump intensity in its dependence on 6 P . All curves 
are normalized by one at their maxima. 

taken in the Gaussian form 
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(||-geomctry) we used the analytical expression of Eq. 
(fT7| . On the other hand, as in the case n' p = the sinc- 
function (Fig. [3ja)) is wider than the pump, the coinci- 
dence curve in the _L-geometry is close to the pump an- 
gular distribution curve in its dependence on (9\ +92)/% 
at a given 6*2- In particular, at 6*2 = the coincidence 
momentum probability distribution in the || -geometry is 
given by 



dw 
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Validity of this equation is pretty well confirmed by the 
calculated coincidence and pump (\E p (9 p )\ 2 ) curves in 
Fig. EJa): the solid- line curve is twice wider than the 
dotted-line one. 

As a whole, for the taken laser and crystal parameters, 
the case n' p = (^-geometry) does not provide conditions 
for observing high entanglement of the arising biphoton 
state. Indeed, the coincidence and single-particle distri- 
butions are seen in Fig. Ela) to be rather close to each 
other. Their widths are equal to A9 < f[ w 8 mrad and 
A9[ s ]_ ks 12 mrad, which corresponds to the width ratio 
R k x = A9\ s l/A9[ c) ± w 1.5. In reality the width A9\ s) 
and the ratio Rk j_ can be slightly higher than estimated 
here because of the non-collinear phase matching pro- 
cesses to be discussed below. But this increase is not 
too high, and the main conclusion remains valid: the _L- 
geometry does not provide conditions for observing high 
entanglement that can be accumulated in the biphoton 
state under consideration. 

This high entanglement can be seen most successfully 
in the || -geometry. The corresponding coincidence and 
single-particle momentum distributions are shown for 
this case in Fig. [3^6). As in the ||-geometry (n' ^ 0) 
the sinc-function in Eq. (jTTJ) is much narrower than the 
pump, the coincidence angular distribution is identical in 
this case to the narrow peak of the sinc 2 -function located 
near 9\ = and described in Fig. Oc): 



dw 



(c). 







(20) 



The width of this peak found above gives us the coinci- 
dence width of the byphoton angular distribution in the 
||-geometry: A9^ = 0.5 mrad. The single-particle angu- 
lar distribution in the | -geometry is determined by Eq. 
(fT7|) and is shown in Fig. [3^6) in the dashed-line curve. 
The width of this curve is equal to A9^ s j = 47.3 mrad. 
The ratio of these widths is the entanglement-parameter 
for measurements in the | -geometry 



with a = 0.004114. 

In the case n' p = (_L-geometry) integration in Eq. 
(|12[) is performed numerically, whereas in the case n' p ^= 



Rh\\ = 



A9 



A9 



(c) 



47.3 

~olT 



94.6 > 1 



(21) 
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In contrast to the earlier considered case of the _L- 
geometry, the || -geometry of measurements does provide 
conditions for observing very high degree of entanglement 
accumulated in the same biphoton state and not seen so 
well in other geometries of measurements. 

Related to the estimate of the entanglement parameter 
Rk || (EH) j the two new specific effects predicted for the 
1 1 -geometry and differing it from the traditionally con- 
sidered _L-gcometry arc: (1) a very strong narrowing of 
the coincidence and (2) broadening of the single-particle 
distribution curves. The narrowing effect is very strong 
indeed: the width Afl^j is 16 times smaller than the 

(c) 

same width A6\ j_ occurring in the perpendicular geom- 
etry, and 8 times smaller than the pump angular width 
a. This shows, in particular, that the coincidence an- 
gular distribution appears to be narrower than that of 
the pump (contrary to traditional expectations based on 
consideration of the _L-gcomctry scheme). Explanation 
of this effect is related to drastic changes in the structure 
of the sinc-function in Eqs. © and (TTTj) (see Fig. [5]), 
which occur when we move from the _L- to || -geometry. 

As for broadening of the single particle angular bipho- 
ton distribution in the | -geometry, it has a rather sim- 
ple qualitative explanation in terms of the non-collincar 
phase matching down-conversion processes. Their origin 
is illustrated by two pictures of Fig. [U The first pic- 
ture [4(a)] is the same projection of the refractive index 
surfaces as in Fig. HJc). But now the pump wave vector 




(a) (b) 



FIG. 4: (a) Diagram of pump and emitted photon wave vectors 
under the conditions of non-collinear phase matching; dashed lines 
are arcs of circles with radii n /2 and centers at the points O and 
B' , wave vectors are in units of u) v /c. (b) Angular distributions 
of photons emitted in the regime of non-collincar phase matching 
(grey) and of the pump (black), a view along the —Oz' direction, M 
is the point of the maximal achievable non-collincar phase match- 
ing. 

k p is taken slightly deflected from the laser axis direction 
Oz' (upon the angle $ p ). As k p (d p ) is shorter than fc p (0), 
its ending point B' cannot be reached by two collincar 
wave vectors k\ and fe, and these vectors must be non- 
collinear both to each other and to k p {d p ). As it's seen 
from Fig. |3Ja), the wave vectors k\ and ki providing ful- 



filment of the exact phase matching condition are emit- 
ted under angles (with respect to the Oz' axis) exceeding 
•d p . This explains why the non-collinear phase matching 
emission broadens the single-particle probability density 
curves. 

By rotating the wave-vector diagram of Fig. BJa) 
around the direction of the pump wave vector k p we get 
cones along which the non-collincar-cmittcd photons can 
propagate. The cone opening angle changes with chang- 
ing -dp by forming in such a way the total emission area. 
Schematically this area of emission (in the angular space) 
is shown in Fig. |4j6) (shaded in grey). This area is seen 
to be wider than that of the pump (shaded in black) and 
it's seen to be wide both in the O'x' and O'y' directions. 
The last case corresponds to detection of photons in the 
_L-plane. Hence, as said above, emission of photons in 
the non-collincar phase matching regime can make the 
single-particle probability density distribution somewhat 
wider than that described by the dashed-linc curve in 
Fig. ma)- 

The point M in Fig. \M,b) indicates the maximal achiev- 
able angle i? p evaluated, e.g., at a half- width of the pump 
angular distribution, d p = a/2. This value of t9 p provides 
the maximal emission cone opening angle, which is char- 
acterized by the circle in Fig. 0J6) bordering the grey 
area and having its center at M . It's seen in Fig. |4j6) 
that the emission are area is asymmetric with respect to 
the point O' (x 1 and z 1 line crossing) in the O'x' direc- 
tion. This agrees with the structure of the single-particle 
probability density distribution in the case of detection 
in the ||-plane [the dashed-line curve of Fig. [3^6)]. The 
reason of asymmetry is evident: deflection of the pump 
wave vector k p from the Oz' axis to the direction opposite 
to that shown in Fig. BJa) makes k p longer than k p (Q). 
For such wave vectors, in contrast to the case considered 
above, the exact phase matching condition cannot be ful- 
filled at all and, hence, such deflections give almost zero 
contribution to the integral in Eq. p7[) determining the 
single-particle probability density distribution. 

Note that the role of photon emission in the regime of 
the non-collinear phase matching can be diminished by 
special profiling the pump angular distribution. In terms 
of Fig. [T] geometries, if the pump angular distribution is 
made narrow in the O'x' direction but, still, remains wide 
enough in the O'y 1 direction, then the point M in Fig. 
0J6) approaches O' and the radius of the grey area de- 
creases, thus making single-particle angular distributions 
narrower. In the extreme limit of the pump angular dis- 
tribution very narrow in the O'x' direction, in the case 
of photon detection in the _L-geometry we return to the 
results described by the curves of Fig. [3ja). Such an 
asymmetric profiling of the pump angular distribution 
can be provided by focusing with a cylindrical lens or by 
a slit. 

So, the used above simplification that the essential 
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pump wave vectors k p belong to the same plane as ki 
and k 2 can be not too good for describing the single- 
particle angular photon distribution in the _L-geometry. 
But this simplification does not affect both coincidence 
distributions at arbitrary orientation of the observation 
plane (z'£) and any distributions in the || -geometry. 



the crystal length L. Its features are fully controlled by 
the pump angular divergence a. 

The biphoton coincidence coordinate distribution 
dw^ / dx[\ x ' 2=0 oc l^coord^ii 0)| 2 calculated with the 
help of Eq. (j2"4")l is plotted in Fig. [5] in its dependence on 



Coordinate representation 

The biphoton wave function in the coordinate rep- 
resentation can be obtained from the momentum- 
representation wave function by means of a double 
Fourier transformation 

mom 

xejq>[i(&fci f + 6afee)] ) ( 22 ) 

where £i and £2 are coordinates along the observation di- 
rection 0'£ for photons (1) and (2) and $ m0 m(fcu, fef) 
is the same function as &(ki£, k 2 ^) of Eq. |(5J). Inte- 
grations in Eq. (|22[) can be carried out partially in a 
rather simple way only in the case of a sufficiently large 
refractive index derivative. So let us consider only the ||- 
geometry when the detection direction 0'£ coincides with 
O'x' . The integration variables k\ x i and k 2x i can be sub- 
stituted by the angular variables 9 p = {k\ x i + k 2x i)/k p 
and 9 = [k\ x ' ~ k 2x i)/k p °\ In these variables the sinc- 
function in |(5J) takes the form convenient for its subse- 
quent approximation by the 5-function: 



sine 



Lh 



(0) 



, n'9„ , 
2n a \ p p ' 2 



oc <5 [9 V + 



2ni 



(23) 



In this approximation and with the pump amplitude E p 
taken in the Gaussian form [Eq. (JT5J)] , Eq. (f2"2"| yields 



coord 



(x'i, x' 2 ) oc / d9 exp 



ln(2)6» 4 
~2a 2 n'J 



4n' 



(24) 



The (^-function approximation of (|23p and the final ex- 
pression of Eq. (f2~4"|) are valid the sine-function on the 
left-hand side of Eq. (|23p is narrower than the pump 



E p (9 p ) of Eq. (fig]), which gives 




FIG. 5: Coincidence probability density in the coordinate repre- 
sentation in its dependence onij kp /2 (solid line) ; dependence of 
the bipartite probability density on (11 + X2)k p °^ /2 at x\ — X2 = 
(dashed line) and on (x\ — X2)k p °" > /2 at xi + £2=0 (dotted line). 



the dimensionless product x[k p /2 (the solid-line 
curve). In the same picture the probability density 
I* 



'(x^x's) is plotted in its dependence on £+ = 
{x[ + x' 2 )k p 0) /2 at = (x[ - x' 2 )k p 0) /2 = (the 
dashcd-linc curve) and on £_ at £ + = (the dotted- 
line curve). The widths of these curves are equal to 
Ax[ {c) k p 0) /2 = 88, A£+ = 356.4, and A£_ = 44. The 
distribution is seen to be much wider in the than in 
£_ direction. 



coord 



We did not calculate (yet) the single-particle prob- 
ability density in the coordinate representation. But 
we can use the found width of the coincidence coor- 
dinate distribution for evaluating the EPR-rclated pa- 
rameter of entanglement. Indeed, by taking into ac- 
count that Aki x > = Af?i x kp /2, we get for the ||- 



geometry: Ak\ 



(c) 



x Ax[ (c) 
' (c) _ 



= 0.0005 x 88 = 0.0044 and 

Cepr = l/Afc^, xAxj = 22.7. By taking into account 
that in the no-entanglement case for Gaussian functions 
and for our definitions of widths as the half- height widths 
C^ R = (41n2)-\ we find finally the the EPR related 
parameter of entanglement for the || -geometry is evalu- 
ated as 



Lki 0) \n' p \ 2 
a 



2n 



(25) 



With L, a, and kp we use in this paper the left- and 
right-hand sides of this inequality are equal to 1.1 x 10 4 
and 486, correspondingly, and the condition of Eq. (|25|) 
is pretty well satisfied. Under the condition (|2"5"|) the co- 
ordinate biphoton wave function (|24[) docs not depend on 



EPR 



(0) 
EPR 



= 22.7 x 4 In 2 = 63. 



(26) 



Though less than 80, this parameter is also rather large 
and indicates clearly that SPDC biphoton states under 
consideration are rather highly entangled. 
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EXPERIMENT 
Experimental setup 

The experimental setup is shown on Fig. 3. To generate 
the entangled photons we use type I and 15 mm-lcngth 
lithium-iodate crystal pumped with a 5 mW cw- helium- 
cadmium laser operating at 325 nm. The correlated pho- 
tons generated via SPDC process with equal polarization 
and wavelength 650 nm are separated from the pump 
by dichroic mirror (DM). Interference filters centered at 
650 nm with a bandwidth of 10 nm are placed in each 
arm of Brown- Twiss scheme. To measure coincidence and 
single-photon distributions in the transverse momenta we 
use the lens with focal length F = 62 cm. Two sin- 
gle photon detectors are positioned in focal plane of the 
lens. Such detector arrangement allows one to measure 
the momentum distribution(s) by scanning one or both 
detectors along with certain direction (see below). In the 
most of cases we fix position of the first detector at the 
maximum of count rate and scan another one to register 
both distributions as a function of detector displacement. 
Since detector moves in the focal plane its position (x) 
relates to the angular mismatch {&) as x = Ft&n9. 




FIG. 6: Experimental setup for measuring single and coincidence 
probability distributions. 



Results and discussion 

The main idea behind performed experiment is to 
check formulas (fT7|) . (|2U|) and (|19p for two geometries 
such that detectors are scanned (*) in the ||-plane con- 
taining optical axis, when 0'£ = O'x' G (xz), and (**) 
in the _L-plane, when 0'£ = 0'y'\\Oy (see Fig. [l}. One 
of the key parameters of the theory described above is 
the pump angular width. Originally our He-Cd laser had 
the angular width equal to 1.5 mrad. To see in experi- 
ment how the angular distrubution of the pump affects 
upon biphoton angular distributions, we have artificially 
and anisotropically broadened the pump distribution in 
angles by installing in front of the crystal a slit. As a 
result, the pump angular distribution remains localized 
along the direction parallel to the slit (at the same level 
of 1.5 mrad as was without a slit) but it spreads in the or- 



thogonal direction. The slit was measured to be 40 mkm 
wide, which corresponds to the pump angular width in 
the direction perpendicular to the slit equal to 4.1 mrad. 
Alternatively we have used in some measurements a lens 
instead of a slit to provide axially symmetric (isotropic) 
broadening of the pump up to the same width of 4.1 
mrad. Comaprison of results of such measurements was 
used for evaluating the role of the pump angular broad- 
ening for efficiency of emission processes arising under 
the non-collinear phase matching conditions. 

As mentioned above, in the experiment the detection 
plane was always horizontal. The slit was installed always 
vertically to provide angular broadening in the horizontal 
direction, whereas the crystal optical axis could be lying 
either in the vertical or horizontal planes. In terms of 
notations and definitions used above and introduced in 
the picture of Fig. Ufa) this means that our experimen- 
tal measurements correspond to one of the following two 
situations: (A) detection and pump angular broadening 
directions are along the O'x' axis, i.e., in the ||-plane or 
in the plane containing laser and optical axes and (B) 
both detection and pump angular broadening directions 
are along the O'y' axis, i.e., in the _L-planc, and these di- 
rections are perpendicular to the plane containing laser 
and optical axes. In the case when lens is used for the 
pump broadening it provides equal angular broadening in 
O'x' and and O'y', independently of the detection direc- 
tion. As in the case B the slit-broadened pump angular 
distribution is narrow in the O'x' direction, we expect 
that under these conditions non-collinear phase match- 
ing processes will be suppressed [see the explanantion in 
Fig. \Mfi)] an d the single-particle angular distribution in 
the direction O'y' will be narrower than in the case of 
lens angular broadening. 

Another key parameter is angular derivative of the 
pump extraordinary refractive index n' p near the exact 
phase matching direction Oz' . The table 1 shows this 
value for different crystals available for producing pho- 
ton pairs. It shows that Lithium Iodate is the best can- 
didate because in this crystal the derivative \n' p \ takes a 
rather large value. Besides, the effect of high entangle- 
ment anisotropy discussed in the paper depends strongly 
on the product Ln' p . So, the second reason why we chose 
this crystal is that the sample of lithium iodate can be 
made rather long. 



Crystal 


phase matching angle <po (deg.) 


n' p 


LBO 


51.47 


-0.0270 


KDP 


54.33 


-0.0395 


BBO 


36.44 


-0.1175 


LUOa 


5.97 


-0.1409 



TABLE I: Angular anisotropy parameter n' p for different non- 
linear crystals 

As it follows from theory developed above, at suffi- 
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cicntly high values of \n' \ the coincidence angular dis- 
tribution of biphotons does not depend at all on the di- 
vergence of the pump [see Eq. PU|) ]. Moreover, being 
determined completely by properties of the crystal and, 
in particular, by anisotropy of the refractive index, the 
width of the coincidence distribution can be done even 
narrower than the pump angular width because. At the 
same time the width of the single-photon distributions 
grows up with the pump broadening. Both narrowing 
of the coincidence and broadening of the single-particle 
angular distributions arc factors resulting in a growing 
degree of entanglement. 

The pictures of Fig. [7] show two sets of an- 
gular distributions received in (a) single-particle and 
(b) coincidences measurements, which are plotted to- 
gether for different geometries, with a slit broadened 
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FIG. 7: Angular distributions of (a) single and (b) coincidences 
counts for two orthogonal geometries (_L) and (||). 

pump beam. These pictures illustrate clearly that in 
the 1 1 -geometry coincidence distribution becomes nar- 
rower whereas the single-particle distribution broad- 
ens in comparison with the _L-geometry. The corre- 
sponding ratios are Afcj^/Afcjj — 11 for coincidences and 

Afcj//Afcj| =0.41 for singles. For comparison, the cor- 



and Ak± / Akl*' w 0.25. The difference between theo- 
retical and experimental results is not too high, though 
quite visible. Probably, it has different origin for coinci- 
dence and single width ratios. In the case of coincidence 
counts, we think that the main reason for a difference 
between the theoretical and experimental width ratios is 
related to some external factors owing to which the exper- 
imentally measured coincidence width is larger than the 
theoretical estimate (0.75 mrad in experiment compared 
to 0.5 mrad in theory). As for the single-particle width, 
probably in experiment the corresponding curve in the 
_L-geometry did experience some broadening arising from 
the non-collincar phase matching emission processes, in 
spite of the missing slit-broadening of the pump angular 
distribution in the O'x' direction. In contrast to this, in 
theory presented above in the _L-gcometry case the non- 
collincar phase matching emission processes were com- 
pletely excluded because only pump photons with wave 
vectors belonging to the observation plane were taken 
into account. Hence, compared to the real experimental 
situation, the given above theoretical derivation of the 
single-particle angular distribution in the _L-geometry ar- 
tificially narrows the corresponding curve, and this is the 
reason for the observed discrepancy between theory and 
experiment. But let us emphasize once again: the dis- 
crepancy is not too high, and the degree of closeness of 
theoretical and experimental results can be considered as 
quite satisfactory. 

Figs. [5] and [5] present the same experimental results 
only differently regrouped, which allows one to evaluate 
the experimentally found degree of entanglement. 

Fig. [5] corresponds to the case when detector is scan- 
ning in the plane perpendicular to the optical axis (_L- 
geometry), and the slit broadens the pump angular dis- 
tribution in the same direction. The width of the singlc- 
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FIG. 8: Experimental results: normalized single-particle and co- 
incidence distributions for ^-geometry. Curves with errors are the 
coincidence. The pump is broadened by a slit. 
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particle distribution is 25 mrad whereas the width of 
the coincidence one is 8.4 mrad, i.e., twice wider than 
the width of the pump [in accordance with Eq. (fTO)) ]. 

Afc M 

The ratio Afce ^ J - — = 3, so the degree of entanglement 
is not very high. Note, that in the above-presented the- 
ory we got for this ratio the twice smaller value, 1.5. 
The reason is the same as discussed above: in theory 
the single-particle angular distribution is artificially nar- 
rowed because all non-collinear phase matching processes 
are completely excluded from the consideration. 

The results occurring in the case when both scanning 
and pump broadening occur in the plane containing the 
crystal optical axis (||-geometry) are shown in Fig. [9] 
Here the widths of single- and coincidence distributions 
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served effect remains unexplained: for the same ±- 
geometry the coincidence width in the case of lens- 
broadening appears to be smaller than in the case slit- 
broadening, Afc ( ( 1 c c ) ns) /Afc ( ( ^ t) =0.0031/0.0084=0.37. Note 
also that in the case of lens broadening the values of 
the entanglement parameter Rk appear to be somewhat 
smaller than in the case of slit broadening in both _L- and 
||-geometries: in the case of lens broadening Rk± = 2.33 
and Rk || = 67. But the main qualitative conclusions 
remain the same as earlier: at chosen values of the pa- 
rameters almost no entanglement can be seen in the _L- 
geometry and a very high entanglement can be and was 
observed in the || -geometry. 
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FIG. 9: Experimental results: normalized single-particle and co- 
incidence distributions for ||-geometry. Curves with errors are the 
coincidence. The pump is broadened by a slit. 

are 60 mrad and 0.75 mrad correspondingly. Their ratio 
is 80, which is much greater than in previous case though 
is somewhat less than the corresponding theoretical esti- 
mate Rk || = 94.6 (|2ip . The difference can be attributed 
only to some external factors affecting experiment. In the 
case of the || -geometry the above-discussed non-collinear 
phase matching processes are completely taken into ac- 
count in the theoretical derivation, and the explanation 
of the theory-experiment discrepancies given above for 
the _L-geometry does not work for the || -geometry. 

Fig. [TO] presents similar graphs but measured for 
the case when the laser pump is broadened with a 
spherical lens. This means that the pump broaden- 
ing affects angular distributions in both directions and 
so, it is impossible to distinguish its contribution to 
spreading in ||- or _L-planes. According to the pic- 
ture of Fig0J&), the single-particle distribution mea- 
sured in _L-plane must be wider when the pump angu- 
lar distribution is broadened by a lens compared to the 
case of broadening by a slit. This expectation is con- 
firmed by measurements which give (for the _L-gcometry) 
Ak { ^ na) /AkH t) =0.075/0.025=3. However, another ob- 
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FIG. 10: Experimental results: normalized single-particle and co- 
incidence distributions for (a) ||- and (b) _L-geomctries. Curves 
with errors are the coincidence. The pump is broadened by a lens. 

In conclusion, we have shown both theoretically and 
experimentally that there are two different schemes of 
observing SPDC biphoton wave packets. In the tra- 
ditional and alternative schemes the detector scanning 
is assumed to be performed in the planes, correspond- 
ingly, perpendicular and parallel to the plane contain- 
ing optical and laser axes. Owing to the anisotropy of 
the crystal refractive index for the extraordinary wave, 
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a structure of the coincidence and single-particle bipho- 
ton wave packets observable in these two schemes are 
significantly different. In the alternative scheme, coinci- 
dence wave packets demonstrate a very strong narrowing 
compared to the traditional scheme, whereas the single- 
particle wave packets broaden. All this results in a very 
high degree of entanglement that can be observed in the 
alternative scheme of observations, and cannot be seen 
in the traditional scheme. The entanglement parameter 
determined as the ratio of the single to coincidence wave 
packet widths appears to be as high as about 100. The 
degree of the coincidence wave packet narrowing is shown 
to be so strong that in the alternative scheme of obser- 
vations the coincidence wave packet appears to be much 
narrower than the angular distribution of the pump. 
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